Temperature of a Hamiltonian system given as the effective temperature 
of a non-equilibrium steady state Langevin thermostat 
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In non-equilibrium steady states (NESS) far from equilibrium, it is known that the Einstein 
relation is violated. Then, the ratio of the diffusion coefficient to the mobility is called an effective 
temperature, and the physical relevance of this effective temperature has been studied in several 
works. Although the physical relevance is not yet completely clear, it has been found that the role 
of an effective temperature in NESS is indeed analogous to that of the temperature in equilibrium 
systems in a number of respects. In this paper, we find further evidence establishing this analogy. 
We employ a non-equilibrium Langevin system as a thermostat for a Hamiltonian system and find 
that the kinetic temperature of this Hamiltonian system is equal to the effective temperature of the 
thermostat. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.70.Ln 



Fluctuation-dissipation relations (FDRs) relate dy- 
namical properties of fluctuations in systems under equi- 
librium conditions to linear transport properties of non- 
equilibrium systems through the detailed-balance condi- 
tion [l[ . Representative examples of FDRs are the Ein- 
stein relation, which relates a diffusion coefficient and a 
mobility, and the Green-Kubo relation, which relates cur- 
rent fluctuations and the corresponding conductivities. 

In recent years, the properties of fluctuations and lin- 
ear responses to perturbations have been investigated for 
non-equilibrium states even outside the linear response 



regime, specifically in steady state systems @, H, Q| and 
in aging systems HUH!- Although, we cannot expect 
FDRs to be generally valid outside the linear response 
regime, there have been several relations proposed and 
investigated recently that represent extensions of FDRs 
to systems far from equilibrium [l^, EL E3 ■ 

In Refs. [U[l(3], the violation of FDRs is studied in the 
case of a non-equilibrium one-dimensional Langevin sys- 
tem in which a Brownian particle is subject to a spatially 
constant driving force / and a periodic potential U(x). 
Explicitly, the system studied there is 



dU(x) 
dx 



+ / + £(*), 



(i) 



where x(t) is the position of the Brownian particle, 
is Gaussian noise, 7 is the friction coefficient, and T is 
the temperature of the environment. (The Boltzmann 
constant is set to unity.) In this model, in the linear 
response regime, the Einstein relation 



(2) 



holds, where D is the diffusion coefficient and /id is the 
differential mobility defined as 

((*(*) -s(0) -v s (f)t) 2 ) 
D(f) = km ±^ 4- S M±±1, ( 3 ) 
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FIG. 1: Schematic depiction of the "moving thermostat." The 
thermostat moves at a speed v relative to the Hamiltonian 
system. 



Here v s (f) is the steady state velocity of the Brownian 
particle and is known as the Stratonovich formula [§, [H[ . 
However, outside the linear response regime, i.e. for large 
/, the above Einstein relation does not hold. In such 
situations, as an extension of the concept of temperature, 
it is natural to define the following quantity: 



/w) 



(5) 



Then, outside the linear response regime, we have O =/= T. 
Thus, the introduction of 9 allows us to define an ex- 
tended Einstein relation that applies to non-equilibrium 
steady states (NESS) far from equilibrium, although this 
leads to the question of the physical significance of 0. 
In Refs. [H, [l(| , in order to elucidate the physical signifi- 
cance of 6, a large-scale description of the system was de- 
rived by applyinga perturbation method to the Fokker- 
Planck equation [9( and by considering a finite time av- 
erage of the Langevin equation |10j. With these treat- 
ments, it was found that O plays the role of a tempera- 
ture in the large-scale description of the non-equilibrium 
Langevin system ([T]), and for this reason, it is referred to 
as an effective temperature. 
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FIG. 2: K(t) obtained using a stationary thermostat as a 
function of time in the cases / = (upper) and / = 10 
(lower). These results were obtained from 5000 samples. 




FIG. 3: K(i) obtained using a moving thermostat as a func- 
tion of time in the cases / = 5 (upper), / = 10 (middle) 
and / — 15 (lower). These results were obtained from 5000 
samples. 



In this paper, we present a study that further estab- 
lishes the role of as an effective temperature for NESS. 
Here, we employ a Langevin system in a NESS as a ther- 
mostat for a Hamiltonian system, and we investigate the 
temperature of the Hamiltonian system established by 
this thermostat. More precisely, we set out to determine 
whether the kinetic temperature of the Hamiltonian sys- 
tem is equal to O in this situation. We find that in fact 
the kinetic temperature is equal to in the case that the 
thermostat moves at the speed v — —v s (f), so that the 
average velocity of the Brownian particle relative to the 
Hamiltonian system is zero. (See the schematic in Fig. [1] 
and Moving thermostat for details.) 

Stationary thermostat. In Fig. [TJ we present a 
schematic depiction of the model we study. In this 
section, we consider a combined system consisting of a 
Hamiltonian system in contact with a Langevin thermo- 
stat in the case that the two are relatively at rest, i.e. 
v = 0. 

The Langevin thermostat consists of N Brownian par- 
ticles (N — 20) which are confined to move along a single 
direction, say the x direction. Because there is no in- 
teraction between the Brownian particles, the statistical 
properties of each are the same as those in the model (JTJ , 
studied in Refs. d, Each Brownian particle is sub- 
ject to a constant driving force / and a periodic potential 
U(xt) = (Uo/T)sm(2nxi/£), where x% represents the po- 
sition of the z-th particle. The size of the thermostat is 
chosen as 20£ (— 10£ < Xi < Wi), and periodic boundary 
conditions are imposed on the Brownian particles. 

The Hamiltonian system we consider is a one- 
dimensional system consisting of a single particle. Each 
Brownian particle in the Langevin system interacts with 
this Hamiltonian particle through the potential U m t — 
e(xi — x H ) 2 /2 for \xi — x H | < r c and Ui nt = other- 
wise, where x H is the position of the Hamiltonian parti- 
cle, and r c is the cut-off length of the interaction. This 
particle is confined to the region —5£< ir H < 5£ by 
wall potentials of the forms C/l(£ H ) = (x H + hl)~ A and 



U R (x H ) = {x H - 5£)-\ 

The time evolution of the i-th Brownian particle is 
described by the one-dimensional Langevin equation 



dUjxj) 



+ j 



dU int {xi - x H ) 



dxi 

2jTS(t-t')6 h3 , 



(6) 



and that of the Hamiltonian particle is described by 

du int (x t - x H ) ac/ L (x H ) du R (x R ) 



N 

E 

i=l 



dx 11 



dx^ 



dx^ 



(7) 



In our numerical simulation, the velocity Verlet method 
was adopted to integrate the equation of motion ([7]) with 
a time step At = 5 x 10~ 5 , and we used the parameter 
values T = 1, 7 = 1, I = 1, Uq = 3, e = 1, m = 1, r c = 
4and < / < 25. As to an initial condition, x H (0) = 0, 
Xi(0) = i/2 - 10 (i = 1,2, ■ ■ ■ , 20) and i; H (0) was chosen 
randomly according to a Gaussian distribution. 

We define the velocity fluctuation of the Hamiltonian 
particle as 



K(t)^ m {(v*(t)i)-(v*(t)) 2 } 



(8) 



In Fig. [21 we plot K(t) as a function of time for / = 
and / = 10 in the case of a stationary thermostat. Next, 
we define the kinetic temperature as 



K = lim K(t). 

t — >oo 



(9) 



In the case / — 0, we find that K = 1.00±0.016, which is 
equal to the temperature of the environment (T = 1). In 
the case / = 10, because the Brownian particles exhibit a 
non-zero average velocity maintained by /, significantly 
more heat flows into the Hamiltonian system than in the 
case / = 0. But in this case, as in the / = case, 
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FIG. 4: The kinetic temperature K obtained using a moving 
thermostat as a function of / (circles). The solid curve de- 
notes the analytical solution of 0(/) (= D(f)/ fld(f)) f° r the 
model (HJ) 0. 



FIG. 5: Schematic depiction of the heat conduction system 
using the effective temperature. 



as t increases, K (t) approaches a constant value, with 
a kind of stationary behavior being established between 
the Hamiltonian system and the thermostat. In this case, 
we find K = 6.22 ± 0.22. (We obtain K by averaging 
K(t) over the intervals t <E [60, 150] and t E [200, 300] for 
/ = 10 and / = 0, respectively.) 

Moving thermostat. Next, we consider the case of a 
moving thermostat. Specifically, we study the situation 
in which the thermostat moves at a constant speed of 
v = —v s (f) relative to the Hamiltonian system, where 
v s (f) is the steady state velocity of the Brownian par- 
ticles @. With such a moving thermostat, the average 
velocity of each Brownian particle measured with respect 
to the spatial coordinate of the Hamiltonian system van- 
ishes. Defining yi = Xi — v s (f)t, we can realize such a sys- 
tem by simply replacing U- ln t(xi — x H ) with U- m t{yi — .t h ) 
in Eqs. © and Q. 

In Fig. [31 we plot K(t) as a function of time in the 
cases / = 5, 10 and 15. Note that in order to obtain 
these results for K(t), we use the analytical solutions of 
v s (f) for the model (|T]) [lij ]. From the data plotted in 
the graphs of Fig. [31 we find K — 1.25 ± 0.05 in the 
case / = 5, K = 1.54 ± 0.03 in the case / = 10, and 
K = 1.68 ± 0.06 in the case / = 15. Comparing the 
middle graph of Fig. [3] with the lower graph of Fig. [21 
both corresponding to the case / = 10, we find that the 
value of K obtained when using the moving thermostat 
differs significantly from that obtained when using the 
stationary thermostat. 

Now, let us compare the above results for K with the 
values of for the thermostat. In Ref. @, from calcula- 
tions of D(f) and it was found that D/^d = 1-24 
in the case / = 5, D/fi^ = 1.52 in the case / = 10, 
and D/nd = 1.67 in the case / = 15 for the model |T]). 
Because T = 1, these values of D//.ia indicate that the 
Einstein relation does not hold for / > 5. Comparing 
these values with the values of K computed presently, 
we find that the relation K = D//id holds when e, which 



represents the strength of the interaction between the 
thermostat and the Hamiltonian system, is sufficiently 
small [Hj]. This implies that the kinetic temperature of 
the Hamiltonian system is equal to the effective temper- 
ature given by ([5]) in the case that the Hamiltonian 
system is in contact with the moving thermostat. In Fig. 
21 we compare K with 6 (= D/na) for various values of 
/. It is seen that the relation 

R = e (io) 

holds, to the precision of the numerical computations. 

Interpretation of Eq. U0\ ). In Ref. 0, a large-scale 
description of the probability density for the model ([T]) 
was derived using a perturbation method, and it was 
found that appears as a temperature in a Fokker- 
Planck equation of the coarse-grained probability den- 
sity. Then, in order to further investigate the physi- 
cal properties of 0, in Ref. [13], a coarse-grained de- 
scription of the motion of a Brownian particle was de- 
rived by computing a finite time average of the Langevin 
equation, rather than analysing the probability density. 
This coarse-grained description is given by the equations 
r(X n+1 - X n )/St = F + E n and (E n E m ) St = 2TQ6 m , n , 
where we have F = Tv s , X n = x(t n ) and t n = ndt 
(n = 0, 1, 2, • • • ), and the time interval St is chosen to 
be sufficiently longer than the characteristic time of the 
system. Here, T and F are uniquely determined as func- 
tions of the parameters that appear in the model ([T]) [l(| • 
Then, using the moving coordinates Y n = X n — v s t ni we 
can describe the large-scale motion of a Brownian parti- 
cle by the equilibrium-form Langevin equation 



In the present investigation, choosing the cut-off length 
of the interaction between the Brownian particles and 
the Hamiltonian particle, r c , to be sufficiently large, we 
considered the change in behavior of the system as we 
increase the number of the Brownian particles that inter- 
act with the Hamiltonian particle. In the case that there 
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FIG. 6: Kj as a function of j in the cases / = 5 (pluses), 
/ = 10 (circles), / = 15 (squares), / = 20 (triangles) and 
/ = 25 (asterisks). 



are many Brownian particles, the Hamiltonian particle 
moves slowly enough that its characteristic time is larger 
than St. Because in this case, when we use the mov- 
ing thermostat, the motion of each Brownian particle is 
described by (fTTj) , we obtain the result (fT0)l . 

Heat conduction. As an application of the moving 
Langevin thermostat, we study the heat conduction sys- 
tem described below (see the schematic depiction in 
Fig. [5]). Here, a one-dimensional Hamiltonian system 
consisting of 10 particles is in contact with two ther- 
mostats, a moving, non-equilibrium thermostat of the 
type described above and an equilibrium thermostat. Al- 
though the temperatures of the environments of both 
thermostats are set to T = 1, it is expected that a non- 
zero heat flux will be observed in the Hamiltonian system 
because O (which differs from T) plays the role of the 
temperature in the moving thermostat. 

Let x^ be the position of the j-th Hamiltonian particle 



(j = 1, • • • 10). In our model, the j-th particle interacts 
only with its neighbors (the j ± 1-th particles) through 
the potential U^ t (xf - xf ±1 ) = {l/2)(xf - xf ±1 ) 2 + 
(10/4)(x^ - xf ±1 ) 4 . Then, only the 1-st Hamiltonian 
particle is in contact with the equilibrium thermostat, 
and only the 10-th Hamiltonian particle is in contact with 
the moving thermostat. 

Defining Kj as the kinetic temperature of the j-th 
Hamiltonian particle, in Fig. [6l we plot Kj. It is seen 
that Kj < Kj + i. This is due to the relation O > T. 
Although > T in our model, it has been reported that 
the case < T can also be realized with an appropriate 
choice of the periodic potential U(xi) This implies 
that we could control the direction of the heat flux by 
altering the functional form of U(xi). 

Conclusion. In this paper, we have investigated the 
use of a Langevin system in a NESS as a thermostat to 
establish the kinetic temperature of a Hamiltonian sys- 
tem. Our main results consist of the relation (fTD")) and 
the data plotted in Fig. [51 both obtained with the use 
of the moving Langevin thermostat. Because the phys- 
ical relevance of effective temperatures in NESS [3, LL(| > 
glassy systems (fj H, 0] and biomolecules [|[ is not yet 
fully clarified, we hope that our study sheds more light 
on it. 
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